arXiv: 1506.06914v2 [quant-ph] 1 Sep 2016 


The Coupled Cluster Method and Entanglement 

Peter Levay^ , Szilvia Nagy^, Janos Pipek^ and Gabor SarosP 

^Department of Theoretical Physics, Institute of Physics, Budapest University of 
Technology and Economics and MTA-BME Condensed Matter Research Group, H-1521 

Budapest, Hungary 

^Szechenyi Istvan University, Department of Telecommunications, H9026 Gyor Egyetem 

ter 1 

^Department of Theoretical Physics, Institute of Physics, Budapest University of 
Technology and Economics, H-1521 Budapest, Hungary 

(19. July 2016) 


Abstract: The Coupled Cluster (CC) and full Cl expansions are studied for three fermions 
with six and seven modes. Surprisingly the CC expansion is tailor made to characterize the 
usual SLOCC entanglement classes. It means that the notion of a SLOCC transformation 
shows up quite naturally as a one relating the CC and Cl expansions, and going from the Cl 
expansion to the CC one is equivalent to obtaining a form for the state where the structure 
of the entanglement classes is transparent. In this picture entanglement is characterized by 
the parameters of the cluster operators describing transitions from occupied states to singles, 
doubles and triples of non occupied ones. Using the CC parametrization of states in the seven 
mode case we give a simple formula for the unique SLOCC invariant Then we consider a 
perturbation problem featuring a state from the unique SLOCC class characterized hy J 7 ^ 0. 
For this state with entanglement generated by doubles we investigate the phenomenon of 
changing the entanglement type due to the perturbing effect of triples. We show that there 
are states with real amplitudes such that their entanglement encoded into configurations of 
clusters of doubles is protected from errors generated by triples. Finally we put forward 
a proposal to use the parameters of the cluster operator describing transitions to doubles 
for entanglement characterization. Compared to the usual SLOCC classes this provides a 
coarse grained approach to fermionic entanglement. 

PACS: 02.40.Dr, 03.65.Ud, 03.65.Ta 
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1 Introduction 

According to the current paradigm quantum entanglement can be regarded as a new and 
efficient resource [I]. Suppose that we have a multipartite quantum system consisting of a 
number of subsystems. Entanglement can then be created between these subsystems via 
applying a set of global operations on the full system by switching on certain interactions. 
These can be effected via using successively a discrete sequence of quantum gates, or via 
applying a continuous evolution operator generated by some Hamiltonian. The particular 
type of entanglement obtained in this way is defined by prescribing a certain set of local 
operations acting on the subsystems that will result in equivalent states of the full system. 

The most common operations used in quantum theory are of course the unitary ones 
resulting in unitary quantum gates and evolutions [I]. These are the ones leaving invariant 
the scalar products inherently connected to the structure of the multipartite space of states. 
Choosing the prescribed set as the local versions of such operations gives rise to the definition 
of entanglement classes as the set of equivalence classes of the space of states of the full 
system under the action of the group of local unitary (LU) transformations. Hence, two 
states are possessing the same entanglement if and only if they can be converted to each 
other via some local unitary transformation[ 2 . 
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However, for practical purposes it turned out that this classification scheme is too 
restrictive [3l 11]. First of all, local unitary classification is providing a fine graining to the 
full set of states, whereas sometimes a coarse graining would be more desirable. Further¬ 
more in quantum information processing manipulations of more general types then unitary 
ones are also allowed. These are physical manipulations (including communication between 
parties applying classical channels) resulting in the conversion of states back and forth with 
(generally different) success probabilities. These operations are called SLOCC operations 
(stochastic local operations and classical communication). They can be represented math¬ 
ematically via the local action of the general linear group of invertible transformations of 
which the set of local unitary ones is merely a subgroup. The SLOCC entanglement classes 
are then defined as the equivalence classes of the set of states of the full system under the 
action of the group of local invertible SLOCC transformations. This new set of equivalence 
classes gives rise to a coarse graining of the space of states of the multipartite system. 

Notice however, that unlike for a local unitary transformation the physical meaning for 
a particular local invertible one is not so obvious. Let us take the example of a multiqubit 
system. A local unitary transformation on a particular qubit is represented by a 2 x 2 
unitary matrix which can be implemented by using physical equipment like Stern-Gerlach 
magnets, certain arrangements modelled by 2 x 2 Hamiltonians like the ones featuring the 
interaction of 1/2 spins with external magnetic fields etc. On the other hand though being 



easy to implement by a natural physical setup. 

The notion of SLOCC classes can also be generalized to systems with indistinguishable 
constituents[Sl [Bj. Since the subsystems in this case are indistinguishable the local invert¬ 
ible transformations representing SLOCC manipulations should be identical. One can then 
ask what kind of physical representations of SLOCC transformations are available in this 
context? The aim of the present paper is to show that in the special case of choosing our 
system with indistinguishable constituents to be a fermionic one the notion of SLOCC trans¬ 
formations is intimately connected to a physically well established method used by quantum 
chemists for a long time. This method is the Coupled Cluster (CC) method]?] [Hj- 

In the CC method the usual expansion of the state vector I'k) of the fermionic system 
given as a linear combination of Slater determinants (Cl expansion) is replaced by a new 
one. The CC expansion is constructed via the action of a set of commuting cluster operators 
= 1,2,... on a special Slater determinant l^ko) singled out by special physical consid¬ 
erations. The single particle states comprising j'l'o) are called the occupied ones. Then the 
cluster operators Ti describe single, double, triple etc. transitions from the occupied single 
particles states to non-occupied ones. Hence in this CC picture the state l^*) is regarded as 
a deviation from the distinguished state |4'o) effected by multiple transitions from occupied 
states to non occupied ones. 

In this paper by employing three-fermion systems with six|Bl lTUllll] and seven [T^ single 
particle states (modes) we would like to draw the readers attention to the fact that SLOCC 
transformations show up quite naturally within the framework of the CC method. In this 
context invertible transformations such as the aforementioned ones of upper triangular form 
are having a clear cut physical meaning. We also demonstrate that going from the Cl expan¬ 
sion to the CC one plays the role of obtaining a form for a fermionic state where the structure 
of the SLOCC entanglement classes is transparent. The reason for restricting our attention 
merely to these special cases stems from the fact that apart from these ones (and the much 
more elaborate cases of three and four fermions with eight modes [131 [14]) the structure of the 
SLOCC classes and its invariants and covariants needed for our considerations is not known. 
Recall however, that precisely these special cases were the ones where the Borland-Dennis 
inequalities (TB] |TB] of quantum chemistry have been established. As is well-known by now 
these studies have given rise to recent developments culminating in the introduction of the 
generalized Pauli principle in connection with the famous N-representability problem |17lll8j 
and the idea of entanglement polytopesfUl [20]. In this spirit we conjecture that though 
illustrated merely in the simplest multipartite cases our results should have a generalization 
valid also for multi-fermionic systems with an arbitrary number of modes. We hope that 
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our observations will pave the way for interesting future applications. 

The organization of this paper is as follows. In Section 2. we consider the case of 
three-fermionic systems with six modes. First by giving an explicit dictionary between the 
corresponding expansion coefficients we relate the Cl and CC expansions. Then we look 
at the structure of the unique SLOCC invariant T) of fourth order serving as a convenient 
measure of entanglement. We observe that there is a dramatic simplification in the structure 
of this invariant in the CC picture. In this case no contribution from the cluster operators 
Ti appears. It turns out that the reason for this is the fact that Ti generates SLOCC 
transformations, hence the corresponding term can be removed from the expansion of 
jT). Next we show that the remaining expansion coefficients related to T 2 and T 3 are 
characterizing the four nontrivial SLOCC entanglement classes in a simple manner. Section 
3. is devoted to a study of the seven mode generalization of the results of Section 2. 
After recalling the structure of the basic invariant of seventh order, the covariants and 
the SLOCC classes we again relate the CC and Cl expansions. In the CC picture we 
hnd a surprisingly simple expression for As a demonstration how the entanglement in 
this case is manifested in the structure of the cluster operators T 2 and we consider a 
perturbation problem. As a first step we choose a maximally entangled state defined by the 
condition J ^ ^ which is generated merely by the expansion coefficients of (doubles). 
Then we modify this state by adding also the contributions coming from T 3 (triples). This 
perturbation due to triples is characterized by four complex numbers. As a next step we 
derive a constraint for the triples to be able to induce a transition to a different SLOCC 
class. It turns out that the constrained set of parameters define a six dimensional complex 
manifold which is known as the deformed conifold, i.e. a deformation of a six dimensional 
cone. Taking the perturbation parameters to be real reveals a new effect. If we change the 
sign of certain parameters of the doubles then the perturbation due to triples cannot induce 
a transition to a different SLOCC class. Hence in this special case the maximal entanglement 
encoded into the parameters of the doubles is protected from the perturbating effect of the 
triples. In section 4. we put forward a proposal to use the parameters of the cluster operators 
describing transitions to doubles for entanglement characterization. Compared to the usual 
SLOCC classes this provides a coarse grained approach to fermionic entanglement. Using 
our model systems familiar from Sections 2. and 3. and the new system of 4 fermions with 8 
modes we offer some mathematical evidence supporting our proposal. Our conclusions and 
some comments are left to Section 5. For the convenience of the reader we put some of the 
technical details to an Appendix. 

2 Three fermions with six modes 

2.1 Relating the Cl and CC expansions 

We consider three fermions with six single particle states or modes. Let U be a six dimen¬ 
sional complex vector space and V* its dual. We regard U ~ C® with {e^}, ^ = I, 2,... 6 the 
canonical basis and {e^} is the dual basis. Elements of V will be called single particle states 
or modes. We tacitly assume that U is a Hilbert space also equipped with a Hermitian inner 
product, but we will not make use of this extra structure. We also introduce the twelve 
dimensional vector space 

V = V®V*. (1) 

An element of V is of the form x = v + a where w is a vector and a is a linear form. In the 
case of fermions according to the method of second quantization to any element x G V one 
can associate a linear operator x acting on the fermionic Fock space J' as follows. We define 
fermionic creation and annihilation operators as 

e ^ = h ^, ^= 1,...6 ( 2 ) 

satisfying the usual fermionic anticommutation relations 

{p^,h^} = (50i, = {n^,hi,} = 0. (3) 
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We define the vacuum state |0) G by the property 


n^|0) = 0. (4) 

Then T is spanned by the basis vectors 

1 < </X 2 < • • • </r/c < 6 . (5) 

In particular for fc = 3 an unnormalized three fermion state with six modes can be written 
in the form 

\'>P) = ( 6 ) 

where ipp^p is totally antisymmetric having 20 independent complex amplitudes. 

The SLOCC group is the one of invertible 6 x 6 matrices with complex entries i.e. 
GL{6,C). It is acting on a state via transforming its amplitudes as 

i^p.p^ SgGL{6,C). ( 7 ) 

Let us now split the single particle states to ones that are occupied and not occupied. 
The ones that are occupied will be labelled as 

i,j,k= 1,2,3 ( 8 ) 

and the ones that are not as 

a,b,c= 1, 2, 3. (9) 

Hence we have a new labelling for the six modes as 

{I,2,3,4,5,6}={1,2,3,T,2,3}. (10) 

We take the Hartee-Fock (HF) state as the Slater determinant 

Itpo) =p'^fp^\0). (II) 

We see that this distinguished element of IF is built up from single particle states that are 
occupied. Now the Coupled Cluster (CC) and full Cl expansions are respectively 

1^) = gfi+f2+f3|^o) (12) 


and 


\ip) = (i + Ci + C'2 + C3)|V^o). 


(13) 


Here 

fi=TJrni, f2 = ^Tab^^p‘^n,p% fg = (14) 

and 


Cl = CaP°‘hi, C 2 = ^Cab^P°‘niP^nj C 3 = Cj^^‘^^p'^nip'^n2P^n3 (15) 

are called the cluster operators. Notice that by construction the set of operators {Ti, T 2 , T 3 } 
called singles, doubles and triples are commuting. 

For the expansion coefficients of the ”T”s and ”C”s it will be useful to adopt the following 
labelling convention. First we relate the expansion of Eq.(IB]) in terms of the 20 komplex 
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numbers komplex numbers of the Cl expansion in terms of the coefficients 

of the ”C”s 

1 = ^123, = (16) 

For notational simplicity we rename the 20 numbers of this 1 + 9 + 9 + 1 split in terms of 
two complex numbers a, /? and two 3x3 matrices B as follows 

a=l, B\ = Ca\ /3 = Ci^l23_ 

Similarly for the CCparameters we introduce yet another 1 + 9 + 9 + 1 split featuring the 
complex numbers rj^ ^ and the 3x3 matrices X, Y 

V=h X\ = Y\=Ta\ i = (18) 

In order to complete the dictionaries between the different expansions of Eg. (151) . (IT^ . and 
m we need to relate the sets (a, A, B, (3) and ( 77 , X, Y, ^). 

In order to do this let us also compare expansions m and (US. Then we get 

Ci=Ti, 4=f2 + if2, 4 = fa + fifa + if 3 (19) 

fi=Ci, f 2 = C 2 -\cl fa = 4 - + icf ( 20 ) 

One can check that 

fif2pW|0) = Tr{XY)p~fp^\Q) (21) 


if3pW|0) = (Detr)pWlO) (22) 

6 

iffpWlO) = (F#)f £a6cl5W|0) (23) 

f3pW|0)=epW|0), (24) 

where Y'^ is the adjoint matrix of Y i.e. 

rr“ = (DetY)I (25) 

with / the 3x3 identity matrix. Using these results by virtue of Ea. dTQll we obtain the 
following dictionary between the CC and Cl pictures 

a = 77 = 1 , B = Y, A = Y^+X, /? = DetU + (X, F) + ^ (26) 


where 


{X,Y) = Ty{XY). (27) 

2.2 The quartic invariant 

It is well-known that we have a quartic combination T>{tl}) of the 20 amplitudes of |f) which 
is invariant under determinant one SLOCC transformations]^, i.e. under the action of the 
group SL{6,C) of the form 

^ SgSL{6 ,C). (28) 
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This invariant gives rise to a natural measure of fermionic entanglement which for embedded 
three-qubit states boils down[9] to the three-tangle describing the phenomenon of entangle¬ 
ment monogamy[2l]. In order to introduce this invariant first we define the covariant [221112] 


— J_F^^PlP2P3P4P5„/, ('oq') 

1/ — ‘Pi'PiP 2 rp3P4P5 

which is transforming under SLOCC as 

K\ ^ (30) 

where S' = S~^'^. Then 

P(y/-) = = \k>^uK\. (31) 

Hence it follows that T’(y/’) is a relative invariant under SLOCC i.e. 

P(V') ^ (DetS)2p(V-) (32) 

and invariant under the SL{6,C) SLOCC subgroup. 

We will need the expression for T’(y/') in the parametrization of Ea. dTTl) . It is given by 
the formula [9] 

Vitjj) = A[k^ - +al)etA + pBetB], 2k = a/3 - (H, B). (33) 


This expression displays the parameters (a, A, B, /3) i.e. the ones of the full Cl expansion of 
lip) of (TT^ . Using Eq. (l26)) let us now give its new expression in terms of the CC expansion 
parameters {ri,Y, X,^). The result is the much simpler expression 

V{iP) = -b 4DetX. (34) 

Hence this expression is not featuring the matrix Y at all! 

Looking back at Ea. dMl) we see that the alternative choice 

a' = 1, B' = 0, = X, /3' = C (35) 

obtained by setting F = 0 in Ea. (l26)) immediately yields Ea. dMll . Hence we suspect that 
the states |V') answering the set (a,H, i?,/3) and \tp') answering the one {a', A', B', pi') are 
related by a special SL{ 6 , C) transformation. 

This is indeed the case. The group SL{ 6 ,C) is having a subgroup of upper triangular 
matrices of the form 

A€Matr(3, C). (36) 

In the language of the set (a, A, B, /3) these transformations are of the form[9] 

a' = a, B' = B + aA, A' = A + 2B x A + aA# (37) 

/3'=/3-b (A,A)-b (B,A“)-baDetA (38) 

where 

2Bx A = (B-bA)‘*-B!*-A‘*. (39) 

Now we can see that with the choice 

A = -F (40) 

from the set of Eq. (IMl) we can obtain the one of Eq. (IMI) . The upshot of these considerations 
is that using a SLOCC transformation corresponding to the matrix of Ea. (l36p with parameter 


I A 
0 I 
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given by (l40l) we can get rid of the Ti term in the CC expansion of Eo. lfT^ . Hence in order 
to reveal the interplay between the CCM and the SLOCC classes it is sufficient to examine 
the structure of the state 

I'i/)') = (41) 

with parameters given by (1351) . 

It is easy to see that this important conclusion is valid also in the most general case (i.e. 
the number of modes can be arbitrary, say N). Indeed, since the operators Ti,T 2 and T 3 
are commuting ones we can write 

l^/i) = (42) 

Now it is easy to check that the operator always generates SLOCC transformations of 
upper triangular form. To cap all this our observation is valid for an arbitrary number of 
(say n) fermions so one can write 

1^-) =e^i|4''), 14-') = € GL{N,C) (43) 

where 

|vI/o)=pV---p"|0). (44) 

Hence again |'I') and |4'') are in the same SLOCC class. In summary: in the CC picture 
SLOCC entanglement is characterized merely by the cluster operators T 2 ,... ,Tn. Then the 
CC expansion featuring only these operators represent deviations from the separable class 
represented by a single Slater determinant i.e. a HF state of the form 

|4'B) = e^i|4'o). (45) 

The state |4 'b) can also be written as a single Slater determinant in a suitable basis. In the 
case of our example of three fermions with six modes we have 

\^b) = = if + Y\r){f + + H^:p")|0). (46) 

2.3 SLOCC classes and CCM 

As described elsewhere [9] given a state \4>) one can introduce a dual state \tp) which is a 


cubic expression of the original amplitudes of \'ijj) as follows 

\^)^ia,A,B,P) (47) 

a = 2aK + 2DetB, i = 2(/3H>* - 2H x A*) - 2 kA. (48) 

^ =-2/3k - 2DetA, B =-2{aA^ - 2A x B^) + 2 kB, (49) 

where k and the cross product are defined in Eas. ll^ and (l39ll . One can also introduce the 
following set of quadratic polynomials 

Qi = api - AB, Q 2 = A" - PB, Q 3 = B^ - aA. (50) 


In fact these 3x3 matrices are just coming from the blocks of the covariant of Ea. (l2^ hence 
K,p = 0 iff Qi = Q 2 = Qa = 0. Then it is known[9l [12] that the SLOCC classes can be 
characterized as follows. 

Let us denote by Aig our space of three fermions with six modes. Clearly A4o — A^C®. 
We make the further definitions 

Ml = { 1 /; e MoiVii/;) = 0} (51) 
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7^2 = {V’ e = 0} 


(52) 


7W3 = {V'e7Wol^V-=0} (53) 

Notice that in the case of the vanishing of the covariant automatically implies 
= 0,'i/j = 0. Now we have[S] four nontrivial C?L(6,C) orbits which we call SLOCC 
entanglement classes 

OgHZ = Mo - Ml = [(pl 23 +^123 ^^T 23 (54) 

Ow = Mi-M2 = [(p^"^+p^^ + p^"^|0))] (55) 

Obisep = M2 — Mo = + P^^^|0))] (56) 

OsEP = Mo = [p^‘^^m (57) 


where the notation []'(/))] refers to the SLOCC orbit (equivalence class) of the unnormalized 
state \'ip) 

The SEP (separable) class is the class consisting of states that can be written as a single 
Slater determinant. The BISEP (biseparable) class corresponds to the set of tripartite 
states that can be transformed via SLOCC to states consisting of the sum of two Slater 
determinants containing a common factor of single particle states. The W and GHZ classes 
form the genuine entanglement classes. The names separable, biseparable, W and GHZ are 
originating from the fact that when considering embedded three-qubit systems living inside 
our three fermionic ones pun] these classes correspond to the separable, biseparable, W and 
GHZ classes known from studies of three-qubit entanglement [3]. 

Since 27 is a SL(6,C) invariant and |'0) and K are covariants one can calculate these 
quantities merely for the transformed representative \’ip') also belonging to the orbit of \'4>). 
Using Ea. (l35l) one obtains 

a' = C, A' = -(X, B' = -2X^, P' = -^2 _ 2DetA (58) 


Q'i=ii, Q'2 = xK Q'o = -X. ( 59 ) 

From these results we see that: 

I. The matrix Y plays no role in the SLOCC classes. 

II. The SEP (single Slater determinant) states are precisely the states characterized 
by the constraints ^ = 0 and X = 0. For these states 27, Idt) and are vanishing 
automatically. 

III. The BISEP states are characterized by the constraints = 0, X'^ = 0 , DetX=0 but 
X ^ 0 . These constraints automatically yield 27 = 0 and |^) = 0, however ^ 0. 

IV. The W states are characterized by the constraint -|- 4DetV = 0, and the nonva¬ 
nishing of at least one of the quantities of Ea. ((58ll . So for instance we are in the W-class if 
we have ^ = 0, V ^ 0, V** ^ 0, DetV = 0. 

V. The GHZ states are characterized by the constraint -I- 4DetV ^ 0. These states 
are dense within the full set of states. 

Recall that ^ and X characterize the contributions coming from triples and doubles 
of cluster operators. From the above results we also see that in order to Hnd out which 
entanglement class a state belongs the derived quantities of X'^ and DetV should also to be 
taken into consideration. Summarizing: interestingly the notion of a SLOCC transformation 
shows up quite naturally as a one relating the CC and Cl expansions. Going from the Cl 
expansion to the CC one seems to play the role of obtaining something like a simplihed form 
for the three fermion state where the structure of the entanglement classes is transparent. 
In order to see how our observations generalize in the next section we turn to a more 
complicated system. 

8 


3 Three fermions with seven single particle states 


3.1 Invariants and covariants 

A three fermion state with seven modes can be written in the following form 

I't) = (60) 

where I,J,K = 1,2,... 7 and I'/jx is a totally antisymmetric complex tensor having 35 
independent components. The SLOCC group is GL(7,C) and the group action is the usual 
one 

SgGL{7,C). (61) 

The basic covariants are[T^ 

(62) 

NiJ = (63) 

= ( 64 ) 

Under SLOCC the transformation properties of these covariants are as follows 

^ YiS)S^rS^ j.S'Y'Y^'Y'K, (65) 

Nij ^ YiS)S'/'S'/Nrj, ( 66 ) 

L" {V)etSfS^rS\,'Y^' (67) 

where S' = (5“^)^. Notice that the 7x7 matrices Njj and are symmetric. 

From these covariants one can form a unique algebraically independent relative invariant 

J is A relative invariant meaning that under SLOCC it picks up a determinant factor 

JY) ^ (DetS')3j^(l') (69) 

hence it is invariant under the SLOCC subgroup SL{7,C). Defining 

= (70) 
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we have the alternative formula 

(J(^))3 = DetS. (71) 

3.2 Relating the CC and Cl expansions for seven modes 


Similarly to the six mode case for a coupled cluster description of this three fermion system 
we split the modes to ones that are occupied labelled by i,j, fc = 1,2,3 and the ones that 
are not occupied by a, 6 , c = 1, 2, 3,4. Now the CC and Cl expansions will be just the same 
form as the ones displayed in Eqs. (HHl-dlSl) with the exception of T 3 and C 3 having the new 
form 

T 3 = -^TabYp‘'nip''n 2 p''n 3 , C 3 = —^CabYp°'nip'‘h 2 jfn 3 . (72) 
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It is rewarding to tackle this case as a deviation from the six mode one. We write 


|vI/) = |V;) + |a;), \u:) = (73) 

and |'(/i) is given by Eq.®. The 6x6 antisymmetric matrix w has 15 independent compo¬ 
nents. We write 


UJ = 




(74) 


where E, F are two 3x3 antisymmetric matrices and D is an arbitrary 3x3 matrix. Hence 
in the Cl and CC pictures we group the 35 amplitudes to three 3x3 matrices, and two 
3x3 antisymmetric ones so in the Cl picture we will have five matrices A, B, D, E, F with 
the latter two being antisymmetric, and two scalars a and /3. Similarly in the CC picture 
the five matrices will be denoted as X, Y, Z, U, V with the latter two antisymmetric, and the 
corresponding scalars are rj and The scalars and the first two matrices in both pictures 
have already been defined by Eqs. (HZl and (1181) . The remaining matrices are defined as 

C/ = E,, (75) 


= Fab. (76) 

T/= T^A2^ = Uab (77) 

where in these expressions we have a,b = 1, 2, 3. 

Now similarly to the six mode case we expect that the form of the unique relative 
invariant 77(4)') of Eg. (1681) will be of much simpler form in the CC than in the Cl picture. 
Indeed, in the next subsection it turns out that unlike in the Cl picture where this invariant 
is featuring all of the amplitudes {a,l3,A,B,D,E,F) in the CC one only the quantities 
(77 = 1, X, Z, U) show up. 

In order to prove this again we have to relate the Cl and CC amplitudes using Eq. m. 
Then a calculation shows that 

a = r]=l, /3 = C + Tr(XY) DetY, B = Y, A = X + Y'^ (78) 


D = Z + VY, E = V, F = U+iZ^Y-Y^Z) + [{X+ Y^)v]. (79) 

Here 

E = = [v],j = (80) 

i.e. the latter notation refers to building up a 3 x 3 antisymmetric matrix from a 3 component 
vector. 


3.3 The seventh order invariant 

The data (a, (3, A, B, D, E, F) described in the previous section is directly related to the 
complex amlitudes ijk- They can be used for calculating the covariants and the invariant 
of Ea. (l68l) . Again the data connected to the operator Ti is not needed since it can be linked 
to a SLOCC transformation of upper triangular form. It means that it is enough to do 
the calculation with V = Y = 0. The result for the calculation of the components of the 
covariants Njj and is given in the Appendix. Putting everything together the final 
result is the surprisingly simple expression 

J(^) = -Det(G)-^Det(C/A + CY^), G = hzX + X^Z^), (81) 

4c Z 
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where the matrix G is just the symmetric part of the 3x3 matrix ZX. Notice that the 
expression is regular for ^ = 0 since for arbitrary 3x3 matrices we have 

Det(^ + B)= DetA + + Tt{AB^) + DetB, (82) 

and DetUX = 0 due to = —U. 

Let us find inside this formula our invariant of Ea. (l34l) familiar from the six mode case! 
First notice that due to Eg. (15(11) and the antisymmetry of U we have uu^ = UK Then 
introducing yet another 3 component vector w and its corresponding antisymmetric matrix 
W = [w] as 


W = ^(ZX - X^Z^), - X^Z^)jk (83) 

we obtain 

J(4f) = w'^Gw - ^u^Hu - i (DetZ(f + 4DetX)) - ^^Tr{UXZ'^^), (84) 


H = ^{Z^X^ + X^^Z). 
Here we also used that 


(85) 


BetiZX) = BetiW + G) = DetG + Tr{W*G) = DetG + w'^Gw. (86) 

From Eg. dST]) it is also clear that if we have G = ZX (i.e. when the matrix ZX is symmetric) 
and [/ = 0 then J^(4') = -BetZ{^^ + 4DetX)/4. 

An equivalent way of formulating this can be given as follows. Define a) = and 

apply the restricion Y = V = 0. Then according to Eq.([79|) then D = Z, E = 0 and F = U 
hence one can verify that 


6 j\^) =0 iff G = ZA, (7 = 0. 


(87) 


Recall that Detw = (Pfw)^ where 


Pf(a;) = 


233 ! 




iio ^ 




( 88 ) 


then in the case when V = U = Y = 0 we have DetZ = —Pfw. Using this one obtains the 
result 




(89) 


It can be proved that this result remains valid also in the case when the restrictions V = 
Y = 0 are removed [I^. Hence u}\ip) = 0 of Ea. (l87ll can be regarded as a necessary condition 
for the separability of the seventh order invariant to a cubic and a quartic one. 


3.4 Entanglement classes 

The SLOCC classification of three-fermion states with seven modes is known in the math¬ 
ematical literature as the classification problem of three-forms in a seven dimesional vector 
space[531 [531 [2i]- These results has recently been introduced to entanglement theory [T5]. 
Here we would like to recall the structure of the SLOCC classes. 

Let us consider the state 

|4'_) = (p323 _^123 _^231 _ ^312 _|_ ^5114 ^224 _)_^334 ^|q^^ 

In this special case due to relations dZU-dZi) we have 

(a, A, B, f3, D, E, F) = ( 7 ?, A, A, e, A, U, U) = (1, -1, 0,0, /, 0,0). (91) 
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(Notice that according to our findings at the end of Section 2.2. one can always achieve that 
V = Y = 0 and in this case by virtue of Eos. (17^ and (17^ the CC and Cl labels are just 
the same.) Using now Eg. ((Ml) with Z = —X = I, and ^ = C/ = 0 we obtain 

J(^-) = 1. (92) 

It can be shown[T2] that this three-fermion state with seven modes belongs to a dense SLOCC 
orbit characterized by the property J"('P_) ^ 0. Notice that from the mathematical point 
of view this orbit is similar to the dense SLOCC orbit of the famous GHZ state of the six 
mode case for which we had 2? ^ 0. Indeed apart from some signs the first four terms of 
Ea. (IM|l correspond to a structure already known from Ea. dMl) . 

In order to make the connection with the three-qubit GHZ-state more explicit we intro¬ 
duce the complex linear combinations 


^1.2.3 ^^1,2.3 


^1.2.3 ^ ^1,2.3 _ .^1,2,3^ E4 = zp4_ 

(93) 

Then one can write 


|0) = i (e123_^^T23^ 


(pl23_pl23_ 

f23_^123^ 

(94) 


Then the state on the right hand side of Ea. dMl) corresponds to the structure |000) -I- |111) 
of the unnormalized three-qubit GHZ state. It can be proved that the correpondence is not 
superficial since the three-qubit SLOCC group can be embedded in a consistent manner to 
the fermionic SLOCC one[51 ITU]. Now our state of Ea. dMl) can be given the following form 

i ^22 ( 95 ) 

The reason for the notation |'L_) corresponds to the fact that according to the left hand 
side of Ea. dMl) and Ea. dMl) the V invariant of the GHZ-part is negative (see also Ea. dMl) 
and the preceeding discussion). 

Now introducing the notation 

= ^|0) (96) 

the SLOCC entanglement classes are given in Table 1. We note that for a full character¬ 
ization of these classes other covariants are also needed[Tl]. However, for our purposes it 
is sufficient to consider merely the quantities Njj and J. Notice also that according to 
Ea. d74l) the 6 x 6 antisymmetric matrix defines an alternating (symplectic) form on the six 
dimensional space of single particle states. Then the operator defines a 

canonical form for a nondegenerate symplectic form. This explains the abbreviation SYMPL 
showing up in the second column of our table. The other abbreviations conform with our 
ones familiar from Eqs. O-dSil). 

3.5 Perturbing the canonical GHZ-like state 

In this section we would like to obtain an insight into the interplay between the structure of 
cluster operators T 2 and T 3 and the entanglement classes of the previous section. In order 
to achieve this as a first step notice that in our description of the SLOCC classes for three- 
fermions with seven modes we managed to produce the ” GHZ-like” state belonging to 
the dense orbit merely via an application of cluster operators containing only doubles. This 
means that when building up this state only the expansion coefficients of T 2 were used (see 
Eas. dTO)) . d771) and dM]) L 

Now we would like to perturb this state by adding to it a new term which is also con¬ 
taining contribution from triples i.e. from the expansion coefficients of the operator T3. To 
this end we choose 

|4>_) = -f lx), lx) = (ep^^^ + ulp234_^ ^2^314 ^3^124^ |Q^_ 
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Name 

Type 

Canonical form of 4' 

Rank Njj{'i') 

JW 

I 

NULL 

0 

0 

0 

II 

SEP 

^123 

0 

0 

III 

BISEP 


0 

0 

IV 

W 

^123 _(_ ^123 _|_ ^123 

0 

0 

V 

GHZ 

i^l23 _|_ ^123 

0 

0 

VI 

SYMPL/NULL 

(U^ -k U22 + ^33)^4 

1 

0 

VII 

SYMPL/SEP 

(All -k -f U33)£:4 + ^123 

1 

0 

VIII 

SYMPL/BISEP 

(Ell -k + £;33)^4 ^1(^23 ^23) 

2 

0 

IX 

SYMPL/W 

(Ell -k + £;33)^4 ^123 ^123 ^123 

4 

0 

X 

SYMPL/GHZ 

(U“ -k + £;33)^4 ^123 ^123 

7 

^0 


Table 1: Entanglement classes of three fermions with seven single particle states. Repre¬ 
sentative states of the SLOCC classes are obtained by acting with 'h on the vacuum as in 
Eg. dMl) . The covariant Nij and invariant J are defined by Eas. (l63l) and (l68l) . 


Using Ea. (l55)) we obtain 

= 1 - i + (uT) 2 + (^2)2 ^ ^ (gg) 

Hence we remain in the dense orbit unless the condition 

= 4. (99) 

holds. Notice that the equation 

-f + {u^f = 0. (100) 

defines a six dimensional conifold[55]. Just as a two-dimensional cone is embedded in real 
three-dimensional space as -I- = 0, a real six dimensional conifold is embedded in 

via Eg. dlOOl) . The conifold is a smooth space apart from a singularity at ^ = 

= 0. It is knownpS] that one way to repair this singularity is given by a process called 
deformation under which Eg. dlOOl) is modified to 

-f + {u'^f + (101) 

where Q is a real deformation parameter. Now we see that Ea. dO^ giving the condition 
needed for leaving the dense SLOCC orbit is the one of the perturbing parameters coming 
from triples defining a deformed conifold with the special deformation parameter Q = 2. 

For real parameters this means that if the 4 parameters corresponding to 

the cluster operators describing triples are belonging to a three dimensional sphere of radius 
2 then the entanglement type is changed. The parametrized family of three fermion states 
in this case leaves the dense orbit. From Table. 1. we see that apart from the entanglement 
class corresponding to the dense SLOCC orbit, we have 9 more classe|i|. As a next step we 
show how this class can be identified. 

In order to do this we calculate the covariant iV/j a symmetric 7x7 matrix and look at 
how its rank is behaving as we reach the singular locus of J7($_). Rather than using the 
matrix Njj for convenience we use the one Bjj = —^Njj of Eq. (TTOI) . For the state [Jt-) this 
new matrix becomes just the identity matrix, i.e. we have Bjj = 6ij. Using the expressions 
given in the Appendix for the perturbed state we get 

/ I -U^I + U) -4«\ 

Bjj{^_)= I 0 . (102) 

V 0 1 ) 

^ Taken together with the trivial class represented by the zero state. 
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By applying a suitable (generally complex) orthogonal matrix we can diagonalize B 


^diag ^ g ^ 


V2Q 


QI QI 0 
il-U U-^I y2u I e S'0(7, C) 


^diag _ 


0 0^ 

0 (1 + iQ)/ 0 

0 0 1 , 


-\/20 


where 


Q = V + (ui)2 + [u^Y + [u^Y 


(103) 

(104) 

(105) 


and recall also our conventions of Eq. (1501) . From Eq. dm we see that we are in accord with 
Eg . (1051) and also see that our matrix B fails to be of full rank if the (1991) i.e. Q = 2 condition 
holds. From Eg. (11041) we see that in the degenerate case the rank of B is reduced from seven 
to four. From Table. 1. we se that this reduction of rank of B indicates a transition from 
entanglement class X. to class IX. The former class is the class of the state |'1'_) of Ea. dOOl) 
and the latter is the SLOCC class of the state obtained from Ea. (l73)) by keeping |w) with 
E = F = 0 and D = I and using for \iY) a state from the IE-class of the Eg. (1551) form. 

As an other example let us consider the state 


l^+) = l^+) + lx) 


(106) 


where 

I4/+) = (pl23 + ^123 ^ ^23T ^ ^312 ^ ^ ^224 ^ 

This state is labelled by the set of parameters 

a, A, B, p, D, E, F) = ( 77 , X, F, Z, E, U) = (1, 1, 0, 1, 0, U), (108) 

i.e. we have taken the negative of the matrix X. As a result of this our matrix B takes the 
following form 

/ _j -i(C/-C/) -iu\ 

Bij{<^>+)= \-^{^I + U) I 0 . (109) 

V 0 1 / 

A diagonalization of Bij{^+) yields the eigenvalues 

An2,3 = —^4,5,6 = -^7 = 1- (110) 

Then according to Eq. dm we have J7($+) = — (l -|- jQ^)- Hence for ^ and U real we 
cannot leave the dense orbit. However, for complex parameters if the constraint = —4 is 
satisfied then the rank of B is again changing from seven to four. In order to check this just 
take = (27,0,0,0). In this case the matrix of Eg. (11091) has a block diaginal 

structure containing a 6 x 6 block of the form —(173 -I-7<ti) 0 1 where tri^a are the usual Pauli 
matrices. However 173 -|- 7 CTi is of rank one, hence the 6 x 6 block is of rank three. Hence 
the matrix H(4>+) is of rank four. Then if the constraint = —4 holds the state |'1>+) will 
again belong to class IX. 

The upshot of these considerations is as follows. Taking the perturbation parameters 
to be real reveals an interesting effect. If we change the sign of the X parameters 
of the doubles (see Eas. dOTI) and (11081) ) then the perturbation due to triples cannot induce a 
transition to a different SLOCC class. Hence in this special case the entanglement (which is 
according to Table 1. is of type X. ) encoded into the parameters of the doubles is protected 
from the perturbating effect of the triples. On the other hand if no sign change occurs then 
if the perturbation parameters are constrained to lie on a three dimensional sphere of radius 
two, then the entanglement type is transformed to class IX. Note that the different behavior 
of the states |'f'_) and under perturbations with real parameters is related to the fact 
that these states are equivalent under complex SLOCC transformations but inequivalent 
under real ones [22] ■ 
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4 The role of doubles in characterizing entanglement 

Let us now address the entanglement aspects of the CC method in a somewhat more general 
setting. For an arbitrary number of fermions (n) and single particle states (N) one can use 
the coupled cluster expansion of Eg. (1431) . As we know from this expansion the contribution 
coming from singles can be removed since this merely amounts to a SLOCC transformation. 
Hence for the characterization of entanglement we are left with the the terms featuring 
the exponentials of cluster operators T 2 , T 3 ,... T^. In this section we would like to put 
forward the proposal to characterize fermionic entanglement by keeping merely one of such 
exponentials namely the contribution coming from doubles: T 2 . Explicitly let us consider 
the approximation to |4>) of Ea. (H51) in the form 

1 $) = f2 = ^Tab^^rn^P% ( 111 ) 

where i,j = 1,2,.. .n and a,b = 1,2,... N — n. Naively on expects that for sufficiently small 
n and N the parametrization based on T 2 is capable of recovering all the SLOCC classes. 
At the same time our parametrization of fermionic entangled states in terms of parameters 
of T 2 is reducing the number of entanglement parameters from (^) to ( 2 ) (^2 ”) • However, 
the price we have to pay is that for larger values of n and N the description of entanglement 
in terms of doubles can only produce a coarse graining which is generally insensitive to the 
fine structure of SLOCC classes. At the same time we must bear in mind that though we are 
losing information on the fine deatails of entanglement but at the same time this approach 
could be appealing from the physical point of view. Indeed, entanglement is a resource and 
it is the physical problem at hand that defines the criteria under which we should classify 
this resource. Since most of the interaction terms used in solid state, molecular and atomic 
physics are based on interaction terms with similar structure to that of T 2 our proposal could 
be a valuable tool. Of course in order to check the viability of our approach investigations 
with realistic systems should be considered. Such investigations we would like to perform in 
future works. For the time being in the following we would merely like to offer some solid 
pieces of mathematical evidence in favour of our proposal. 

Let us see first how the usual SLOCC classes are reproduced for our cases of n = 3, 
N = 6,7 in terms of doubles. 

In the N = 6 case according to Eq. ([T 8 |) keeping merely the contribution from the doubles 
amounts to using the 3x3 matrix X for entanglement characterization. According to our 
discussion of SLOCC classes at the end of Section 2.3 we see that in terms of the matrix 
X separable states are the trivial ones with X = X^ = DetA = 0. For biseparable states 
we have X 0 but X^ = DetA = 0. For the classes containing genuine entanglement the 
W-class has A ^ 0 and A^* ^ 0 but DetA = 0, on the other hand for the GHZ-class none 
of these three quantities are zero. This is in accord with Ea. dMll showing that for the GHZ 
class X>($) = DetA ^ 0. Hence in this special case we see that a parametrization in terms 
of doubles intersects all of the usual SLOCC classes. 

In the N = 7 case describing entanglement in terms of doubles means a characterization 
in terms of only the two 3x3 matrices 

x\ = Zka = (112) 

In this case J{^) = —DetG where G is the symmetric part of the matrix ZX. Hence we are 
in class X. of Table I. provided this invariant is nonzero. Glearly the classes H.-V. are the 
ones familiar from the N = 6 case. In this case Z = 0, and the constraints are just the ones 
of the previous paragraph. In order to identify the constraints in terms of A and Z for the 
remaining classes of Table I. one has to look at the ranks of the covariants. In particular for 
the calculation of the rank of the covariant Njj we have to use the expressions as displayed 
in the Appendix where now ^ = 0 and {7 = 0. Notice however, that in order to separate 
classes VI. and VH. further covariants are needed. For the structure of these covariants see 
Ref. [ 12 ]. 

As a less trivial example let us consider some aspects of entanglement characterization in 
terms of doubles for the n = 4, A = 8 case. In this case we have to use the cluster operator 
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T 2 of Eg. dllll) with i,j = 1, 2, 3,4 and a, 6 = 1, 2, 3,4. In order to calculate 1$) of Ea. dllll) 
we calculate 

|4>) = (l + f2 + if|)pl234|Q^ (113) 

with 

T’2|^o) = (114) 


ir||4-o) = \{Tab^^Tj^ - Tab^V^ + 

Since the structure of the entangelment classes in this case is quite involved [13] we are 
content with showing that the orbit which is closed under the SLOCC subgroup SL{8,C) 
(this is a ”GHZ-hke” class) can be recovered by this parametrization in terms of doubles. 
Since almost all orbits are closed with respect to this group action[T3] hence we can conclude 
that for all practical purposes for obtaining a coarse grained picture of entanglement this 
parametrization will do. (It would be interesting to check whether we are still capable of 
discriminating all the fine structure of SLOCC orbits, as we were in the n = 3, TV = 3 case, 
a problem we are not addressing here.) 

In order to show this the only thing we have to notice is that our parametrization with 
doubles obviously intersects the seven dimensional subspace spanned by the vectors obtained 
by acting with the following set of operators on the vacuum. 


Pi =p 


1234 


= 1234 


P2=P 


1313 


=2424 


P3=P 


1212 


-3434 


(116) 


P, = pl324^^2413^ ^ ^1414^^2323^ ^ _pl423_^2314^ p^ ^ _pl234_^3412^ (^^ 7 ) 

According to Ref. m the closed orbits are precisely those that meet this subspace generated 
by the vectors above. Indeed, choosing 


p_12 _ p_34 _ 

^12 ~ ^34 “ 


rj i 34 _ r p 12 _ 7 

^12 — ^34 ~ 


= c, = d, (118) 


?T4 = % = e, = /, + /" = 1 

shows that |$) will be a state in the subspace spanned by P/|0), / = 1,2,... 7. 


(119) 


5 Conclusions 

In this paper by studying simple fermionic systems we established a connection between 
the Coupled Cluster (CC) Method and quantum entanglement. In the CC method for n- 
fermions with N modes one starts with a single Slater determinant comprising a special set 
of n single particle states, called the occupied ones. Then to this distinguished state one 
applies a commuting set of cluster operators of the form e^',i = 1, 2,... n. We have initiated 
a study of the problem how these cluster operators are encoding the SLOCC entanglement 
properties of the system via complex expansion coefficients of operator combinations called 
singles, doubles, triples etc. These operators are describing transitions from the occupied 
modes to the non-occupied ones. Our starting point was the simple observation that the 
operator featuring singles is a special fermionic SLOCC transformation. Indeed, the 
transformation jd^) i—>■ e^^\^) can be rewritten in the form jd^) i—>■ S ® S ® ® S\'^) 
meaning 

SgGL{N,C) (120) 

with S an N X N matrix having the (|551) structure. As a result of this in order to study 
SLOCC entanglement we merely have to look at the structure of the cluster operators 

^2 5 ^3,... Tk. 
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In order to make some progress in clarifying how multipartite entanglement manifests 
itself in these operators we conducted a study on the simplest nontrivial cases of n = 3 
and N = 6,7. In these cases the structure of the SLOCC classes and the structure of the 
invariants and covariants is explicitely known. Surprisingly in the special case of six modes 
the CC expansion is tailor made to characterize the usual SLOCC entanglement classes. It 
means that the notion of a SLOCC transformation shows up quite naturally as a one relating 
the CC and Cl expansions , and going from the Cl expansion to the CC one is equivalent 
to obtaining a canonical form for the state where the structure of the entanglement classes 
is transparent. Using the CC parametrization of states lEas. (fT6ll - (IT^ and Eas. (l75|) - (l^ l 
in the six and seven mode cases we have given simple expressions for the unique SLOCC 
invariants V and J (see Eas. (l34l) and dm) serving as natural measures of multipartite 
entanglement O [Tl] [12]. In the six mode case we have given a full characteriation of the 
SLOCC entanglment classes in terms of the cluster operators corresponding to doubles and 
triples. In the seven mode case we have considered a perturbation problem featuring a state 
from the unique SLOCC class characterized by ^0. For this state with entanglement 
generated by doubles we investigated the phenomenon of changing the entanglement type 
due to the perturbing effect of triples. For the representative of the dense orbit (class X.) 
we have chosen a ”GHZ-like” state generated by doubles. We then used the four complex 
parameters characterizing the triples as perturbations. We have shown that when these 
parameters define a deformed conifold with a deformation parameter Q = 2 then a transition 
to a new entanglement class (class IX.) is possible. We have also demonstrated that there 
are states with real amplitudes such that their entanglement encoded into configurations of 
clusters of doubles is protected from errors generated by triples. 

The results of this paper were primary aimed at drawing the readers attention to the 
existence of the connection between multipartite fermionic entanglement and the structure 
of coupled cluster operators. A natural question to be asked is how our results illustrated 
merely in the simplest nontrivial multipartite (i.e. tripartite) cases generalize to more re¬ 
alistic ones which are containing genuine multipartite systems with an arbitrary number of 
modes. A promising idea in this respect seems to consider entanglement characterization 
via doubles meaning by the entanglement parameters showing up in the cluster operator 
T 2 . Dealing with merely doubles is reducing the entanglement parameters that we have 
to account for, but at the same time we are loosing some of the finer details provided by 
the usual SLOCC classification. However, entanglement is a resource and it is the physical 
problem at hand that defines the criteria under which we should classify this resource. Since 
interactions based on exchange processes involving pairs of modes (particles) are common 
in physics, an appreciation of this idea could provide a basis for further elaborations of our 
proposal. 

In order to fully explore the connection found here between entnaglement and the CC- 
method there are many possible avenues to follow. Obviously the most important of such 
vistas would be a demonstration how the entanglement encoded in the cluster operators can 
be used for studying the properties of realistic (e.g. three electron) systems. Such ideas 
we are intending to explore in future work. Another mainly theoretical line of development 
would be to invoke further results from the theory of spinors [27l US] |29l|3^. In particular 
in order to investigate entanglement patterns one has to use the natural invariants defined 
for spinors [28l|29l|30l|3Tl|T4], and possibly other notions such as the nullity of a spinor [30]. 
Hopefully these structures will provide a further insight into the coarse grained structure of 
entanglement types based on the CC method. 

6 Appendix 

Here for the convenience of the reader for the seven single particle case we give the explicit 
form of the covariants Njj and in terms of the relevant coupled cluster parameters 
(^, X, Z, U). In the expressions below X and Z and U are complex 3x3 matrices where U is 
antisymmetric hence it can also be parametrized as Uab = SabcU^ with a three component 
vector. Recall the index structures Zia, , {Z^)°''‘ and {X^Y^ which are being in accord 
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with expressions like Zia{Z‘^)°‘^ = (I)etZ)Si^ etc. 

iVy = 3[ZX + {ZXf],j, Nab = -3[Z^X* + {Z^X^ f]ab, Nrr = -6DetZ (121) 


Na^ = 3iUX + ^Z^)ar, N„ = 3{Zu),, Na7 = -3e abc{X . (122) 

= -6[(ZX)# + [ZXf^] - 3e^^^e='^-[{ZX)u^{ZX)J^ + {ZX)l^{ZX)iJ, (123) 
L“''= 12[(XZ^“) + (XZ'^“)'^]“''+ 6M“M^ = 6[C^+ 4DetX]. (124) 

L“ = X^-j{Zu)k + 6e*^'=u“(ZX)jfe - 12^(Z‘*)“ (125) 

L*7 ^ Q^£^J>^{zX)jk - 6e^^'^UabX'^jX\, L“^ = -6^^“ + 12e“'’‘=(Z'^X“)6e. (126) 
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